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INVESTIGATION OF NEW RADAR-DATA-REDUCTION TECHNIQUES 
USED TO DETERMINE DRAG CHARACTERISTICS OF 
A FREE -FLIGHT VEHICLE 

By Gerard E. Woodbury and John W. Wallace 
Langley Research Center 

SUMMARY 

An investigation was conducted of new techniques used to determine the complete 
transonic drag characteristics of a series of free-flight drop-test models from principally 
radar tracking data. In addition to fully using the capabilities of both the tracking radars 
and the meteorological measurement systems, preflight trajectory design, exact kine- 
matic equations, and visual-analytical filtering techniques were employed in the study. 

The results of this endeavor were compared with the results obtained from analysis of 
the onboard, accelerometer and pressure sensor data of the only drop-test model that 
was instrumented. 

The accelerometer drag curve was approximated by the radar-data-alone drag 
curve. However, white noise in the radar data prevents a precise definition of the drag 
rise. As a backup method, on the other hand, the radar-data-alone approach has con- 
siderable merit for estimating drag rise, particularly if subsonic wind-tunnel data are 
coupled with the flight data. 


INTRODUCTION 

A series of free-flight drop tests was conducted by the NASA Langley Research 
Center to determine the transonic drag characteristics of a low -drag supercritical body 
of revolution. The drop-test series consisted of three uninstrumented models and one 
instrumented model. The instrumented model contained longitudinal accelerometers and 
various pressure sensors. These instruments were the primary source of data for the 
determination of the drag of the test configuration. The free-flight drag characteristics 
obtained from analysis of these data are presented in reference 1. 

In the planning stages of the test program, meanwhile, it was felt that while the 
instrumented model was being readied, three uninstrumented models could be quickly 
constructed, aircraft dropped, and tracked by radars. By using the radar-data reduction 
methods of the past, these uninstrumented drop tests could provide a quick look at the 



free-flight subsonic and supersonic levels of drag. The tests could also help perfect the 
aircraft-drop and radar -tracking techniques; furthermore, in the event of a malfunction 
of the instrumented model, they could also possibly serve as a backup for the determina- 
tion of the drag rise. It was later decided to explore this latter possibility more thor- 
oughly when, as a result of a stringent accuracy requirement placed on the test series, 
much emphasis was put on using the full capabilities of the tracking radars and the 
meteorological measurement systems. This emphasis coincidentally increased the 
chances of successfully determining the complete transonic drag characteristics of the 
model from radar data alone. 

The purpose of this report is to present the new techniques devised for this effort 
and the results of the application of these techniques to the radar data of the drop-test 
series. The approach taken in this work involved; (1) updating the equations for deter- 
mining the drag from good approximations to more exact expressions written in a form 
conducive to accurate smoothing of the data, (2) using a smoothing technique new to this 
field of endeavor and using it in conjunction with a visual display in order to assure phys 
ical relevance, (3) designing the flight trajectory to maximize the contribution of data 
known to be the most accurate, and (4) employing the most recently developed meteoro- 
logical measurement systems. 


SYMBOLS 


Cp drag force coefficient, — g 

Crp total aerodynamic force coefficient; for small angles of attack C T ~ C D * C X 

C x axial-force coefficient 

d^ base diameter of flight vehicle, cm 

d maximum diameter of flight vehicle, cm 

max 

D drag force, N 

F^(0)j total aerodynamic force relative to Earth-centered inertial axis system, 

N (see ref. 2) 

2 

g sea-level acceleration of gravity, m/sec 

h geodetic height of flight vehicle, m 


2 



M 

< 1=0 

R 

S 

t 

V c 

V E 

V w 

v« 

w 

x 

y 

y E, P 

r E,y 

V 

V 

p 

p 


mass of flight vehicle, kg 
Mach number 

2 

free-stream dynamic pressure, N/m 
flat-Earth horizontal range, m 
reference area, m^ 
elapsed flight time, sec 
atmospheric speed of sound, m/sec 

magnitude of Earth-relative velocity of flight vehicle, m/sec 
magnitude of wind velocity, m/sec 

magnitude of free-stream velocity of flight vehicle, m/sec. 
weight of flight vehicle, mg, N 
north-south component of R, m 
east -west component of R, m 

Earth-relative flight-path angle in vertical plane, deg 
Earth-relative flight-path angle in horizontal plane, deg 
total angle of attack, deg 

elevation angle measured by radar above its horizontal plane, deg 
slant range measured by radar in line of sight to flight vehicle, m 
rate of change of slant range, m/ sec 



p D p measured directly by the radar operating in Doppler mode, m/sec 

2 

O free-stream density of atmosphere, kg/m 

• oO 

xjs azimuth angle measured from north by radar in its horizontal plane, deg 

\h wind azimuth, referenced to north, deg 

A dot over a symbol denotes differentiation with respect to time; a double dot over 
a symbol denotes a second differentiation with respect to time. 

ANALYTICAL METHODS 
General 

A review of the method typically used in the past for obtaining C D from radar 
data alone was made. It was found that this technique had been devised in the days before 
high-speed computers and was initially set up to be done primarily by hand and had nec- 
essarily used many simplifying assumptions. These assumptions are presented in 
appendix A together with the equations based on them for determining C D . Although 
most of the handwork was in time eliminated from the past method by one means or 
another, its equations changed little or not at all. Because of a stringent accuracy 
requirement, preflight calculations were made to check the accuracy of the equations 
when applied to this particular test series. Using the drag-force equation (Al) with a 
computer-generated nominal trajectory of the test configuration, a 3-percent error in 
C D was found in the test region of interest. Thus, it was deemed necessary to devise 
new formulations which are physically more exact than the old ones and to make use of 
the latest technology and equipment. The new radar-data-reduction technique presented 
herein employs a computerized curve-fitting program new to this work in its filtering 
process and uses it in conjunction with a control console. This arrangement permits, in 
all phases of the filtering, visual inspection of the smoothed and differentiated data; 
through direct access to the computer program, it also provides a more expeditious 
means of managing the data. 

The equations for computing position, velocity, and acceleration are written in an 
oblate -spheroidal Earth -centered inertial frame. In addition to being more exact, these 
equations were specifically tailored so that the basic radar measurements (p,i^ r ,0 r ,Pp), 
their time rate of change, and where applicable their second time rate of change would 
be input directly to determine C D . This approach enables due regard to be given to 
the relative accuracy of each measurement by assuring that very accurate data would 
not be degraded by less accurate data in the filtering process, 
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Hand -Calculation Method 


For convenience, the past method may also be referred to as the hand -calculation 
method. No one reference could be found which satisfactorily gave a complete descrip- 
tion of its original formulation. A version of it was reconstructed, however, and is pre- 
sented in appendix A along with its basic assumptions and equations. Basically, two 
radars were involved in the test setup. A continuous -wave Doppler radar positioned 
within a short distance of the launching site measured p^, the rate of change of slant 
range to the model. A remote second radar of a positional nature provided the spa- 
tial location of the model in spherical coordinates ^o,i// r ,0 r ^. (Today's FPQ-6 radar 
combines both of these modes of tracking in one facility.) In general, there were two 
methods of launching the flight models: by rocket propulsion and by a helium gun. As 
the latter name implies, small models were accelerated to low supersonic speed by a 
controlled expansion of the gas (helium) acting on a push-plate within a 6-inch rifle bar- 
rel. The determination of the transonic drag variation was made for both launch modes 
in coasting flights as the models decelerated from low supersonic to high subsonic 
speeds. 

The hand -calculation method is still applicable today. It yields its best results 
when applied to flight tests wherein large drag forces overshadow those forces neglected 
by the simplifying assumptions of the method. It was the neglect of the Coriolis force 
which, in fact, led to the aforementioned 3-percent error in of the low drag bodies 

of this test series. The hand -calculation method could be particularly useful to those 
who are without sophisticated software and who are not laboring under the constraint of 
extreme accuracy as herein was the case. The hand-calculation method was updated 
from time to time primarily to eliminate its handwork and to adapt it to other modes of 
free-flight testing. The results of its application over the years to many flight test mod- 
els are summarized in reference 3. 

New Filtering and Data-Reduction Method 

A block diagram summarizing the new procedures employed in the analysis of the 
flight data is shown in figure 1. The data are first smoothed (filtered) by a cubic spline 
approximation computer program (ref. 4) utilizing a cathode -ray tube console (CSA/CRT). 
On the console, the data are visually presented and checked for wild points, dropouts, 
and any other data anomalies. The data are then fitted to cubic equations (splines) 
between finite boundaries (knots). The degree to which the fit approximates the data is 
also checked visually on the console. The obvious advantage of this technique is that 
engineering judgment based on experience may be exercised at this point, and if not 
satisfied, iterations may be performed varying the number of splines and location of 
knots. The fitting equations are constrained such that the function, its first derivative, 
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COEFFICIENTS OF FITTED 
CUBIC EQUATIONS 



Figure 1.- Block diagram illustrating the flow of raw data to reduced coefficients. 
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and second derivative must be continuous at the boundaries. The first and second deriv- 
atives of the function (data) fitted by cubic equations, naturally, take the form of quad- 
ratic and linear equations, respectively. Thus, it was not surprising to find that while 
the derivatives were continuous at each knot, the derivative time histories were not nec- 
essarily smooth in appearance overall. This was particularly true of the second deriva- 
tives. For that reason, the internally computed first derivatives were also fitted and 
new first and second derivatives were obtained which were more accurate and less erratic. 
The effect of this successive smoothing procedure on the data and its initial first and 
second derivatives is illustrated more thoroughly later. 

All of the equations for determining from radar data alone were incorporated 

into a second computer program which is hereinafter referred to as RADFAST and are 
presented in appendix B. The output of the CSA/CRT program takes the form of coeffi- 
cients of the cubic equations. This output then serves as input to RADFAST where the 
coefficients are used to generate time histories of the nine quantities (p, i// r , and 

their first and second derivatives) that are required by the kinematic equations of 
RADFAST in order to determine displacements, velocities, and accelerations in an iner- 
tial axis system set in the center of a rotating oblate spheroid. The important output 
from the computations of the kinematic equations (eqs. (Bl) to (B38) in appendix B) are: 

(1) Geodetic altitude (eq. (B28)) is required to relate environmental parameters to 
the flight trajectory. 


(2) Free-stream velocity (eq. (B36)) is obtained by adding the wind velocity vectori- 
ally to the Earth- relative velocity. 


(3) Total inertial acceleration (eq. (B19)) is required in the determination of the 
total aerodynamic force acting on the body. 


The total aerodynamic force of the body is computed by the product of the mass of 
the body and the total aerodynamic acceleration in equation (B43) wherein the gravita- 
tional acceleration of the body is subtracted vectorially from the total inertial accelera- 
tion determined from radar data. The total aerodynamic force coefficient of the body is 
computed by 


Crp = 


Ff(Q)j 


q«,s 


1 9 

where q^ = ± p V* . The results of preflight studies presented in appendix C indicated 
that the model could be expected to fly at very small angles of attack (rj * 0) during the 
data period. Thus 
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Finally, with M found by 



the variation of with M may be determined. 

APPLICATION OF TECHNIQUES TO FLIGHT TESTS 

The flight test series was conducted at the NASA Wallops Flight Center. The 
CSA-RADFAST techniques described in the previous section were applied to the radar 
flight data of three of the models of the drop-test series. Although each of the three unin- 
strumented models (designated Rl, R2, and R3) and the instrumented model (R4) were 
successfully dropped and radar tracked, the FPQ-6 radar was not able to obtain any 
Doppler data during the flight of the second model, R2. Consequently, because of the 
incompleteness of the data, the CSA-RADFAST techniques were not applied to the R2 
flight. Since the tests were largely repetitious in nature, the atmospheric and radar 
flight data of only one model are presented herein and may be thought of as being typical 
of that of the three models analyzed. 


Test Models 

The four test models were bodies of revolution with identical exterior contours and 
were designed to have a drag divergence near Mach 1. A sketch of the model is shown in 
figure 2. A photograph of the model mounted on the drop aircraft is shown in figure 3. 
The model length was 114.3 cm and the maximum diameter was 12.7 cm. The fins were 
swept 45° and had biconvex cross section with a thickness-to- chord ratio of 0.03. The 
models were mass balanced to zero, within the tolerance of the balance machine. Each 
model had a mass of approximately 36.4 kg and contained a beacon to aid in radar acqui- 
sition and tracking. 


Test Description 

In each flight test the model was dropped from a T-33 aircraft flying horizontally 
with a nominal altitude of 9144 m and Earth-relative velocity of approximately 183 m/sec. 
The trajectory was directed at the Wallops FPQ-6 radar such that pj-j, the only velocity 
component measured directly and, hence, the most accurate component of the total veloc- 
ity, was maximized. For the same reason the line of sight from the radar was designed 
to be tangent to the flight path at a point in time when M = 1. This design condition was 
not met in the actual flight tests because of Wallops range safety limits. For example, 
in the actual flight test of R4 at M = 1 the flight-path angle y E was approximately 
-61° while the radar elevation look angle 9 r was approximately 20°. 
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Figure 3.- Photograph of flight vehicle on the launch aircraft. 
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Atmospheric Properties 


Measurements of free-stream temperature, pressure, and relative humidity were 
made by rawinsondes for each flight. Free-stream density and speed of sound were 
derived from these measurements. The measured and subsequently derived atmospheric 
quantities of the R4 flight are presented in figure 4. The winds were measured by the 
Jimsphere technique (ref. 5). Briefly, the Jimsphere technique is a method developed 
for measuring winds which employs precise radar tracking of a smooth, superpressure, 
spherical balloon modified by the addition of conical surface-roughness elements. The 
winds are averaged over intervals of 25 m as opposed to approximately 300 m as is done 
in the commonly used radiosonde system. 

Test Data 

The radar data ^p,^ r ,6 r ,pj^ measured during the flight of R4 is presented in 
figure 5. Superimposed on the data is the smoothing of data by the CSA/CRT technique. 
Using the coefficients of the cubic equations fitted to the data in a least-square sense, 
the first and second derivatives of the xp T and 0 r data were generated and are pre- 
sented in figures 6 and 7, respectively. The <9 f plots are particularly illustrative 
of the need for the aforementioned successive smoothing. Its second derivative Q r 
clearly shows in figure 7 that while the derivatives are continuous at each knot, their 
overall time history is oscillatory in nature, apparently about some mean. The dashed 
curves (figs. 6 and 7) represent the new first and second derivatives obtained by fitting 
new sets of cubic equations to the initial internally generated first derivatives of the 
and data. Notice that this procedure does not appreciably change the general value 
of the original first derivatives but does considerably alter the time history of the second 
derivative to a more plausible curve. The successive smoothing procedure was only 
required for the xf/^ and 8 r data since the first derivative of p is measured directly 
by the radar and is the so-called Doppler velocity pp. 

RESULTS OF APPLICATION 

The results of applying the CSA-RADFAST method to the radar flight data' of the 
test models were obtained and then compared with the results of the analysis (ref. 1) of 
the accelerometer and pressure sensor data of the instrumented flight model, R4. The 
two sets of R4 results are of particular interest, for here using different types of data 
from the same flight (a flight in which the model obviously flew the same trajectory, 
encountered the same winds, etc.) it was possible to make a one-to-one comparison of 
their resultant drag variations. 


n 



Terop erasure, K 





(c) Density. 



(e) Relative humidity. 


Pressure, kN/m 


(b) Pressure. 



Speed of sound, m/sec 

(d) Speed oi sound. 


Figure 4.- Meteorological data. 
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Figure 7.- Generated first and second derivatives of the elevation data. 



CSA-RADFAST Method Results 

A total axial-force coefficient was obtained from reference 1 by algebraically 
adding the base and forebody axial-force coefficients presented therein. Then by noting 
that for small angles of attack = C^, the drag variation used herein for comparison 
purposes was realized. 

The results of computing the drag coefficient versus Mach number using the pre- 
viously described techniques are shown in figure 8 and are compared with the results of 
reference 1. It can be seen that the results of the R1 and R3 flights are inconclusive. 
They each have a drag rise at about the same Mach number as that of reference 1, but 
neither approximates the drag curve of reference 1 in appearance. On the other hand, 
the drag curve of R4 does approximate that of reference 1 in appearance. A small ampli- 
tude perturbation on the R4 curve, however, makes it somewhat difficult to precisely 
define its drag' rise. It should be noted here that in this updated method no smoothing was 
performed on anything but the basic radar data. The resultant R4 drag curve could have 
been further smoothed using a priori knowledge of the wind-tunnel drag variation. Having 
a high level of confidence in the subsonic portion of the wind-tunnel drag variation, the 
perturbation in the subsonic radar drag variation would have been eliminated and the 
drag rise could be easily estimated. Thus, as a backup method for determining the drag 
rise, the radar -data- alone approach has considerable merit. 

Additional Comparison 

Independently of the techniques reported herein, an investigation was made by the 
NASA Wallops Flight Center of a Kalman Filter computer program (KAPPA) applied to 
radar tracking data for the purpose of estimating errors in aircraft and missile trajec- 
tories (ref. 6). One set of tracking data used in that study was, again, that of R4. In the 
reference, R4 is designated Supercritical Body No. 3. The only result presented therein 
is a plot showing a variation of drag acceleration with Earth -relative velocity. In order 
to be able to time correlate that result with the CSA-RADFAST results, time histories of 
trajectory parameters were obtained from Wallops via private communication. A com- 
parison between the trajectory results of the two programs is shown in figure 9. It can 
be seen that there is very good agreement between the CSA-RADFAST and KAPPA posi- 
tion and velocity results. Since the model is one and the same, it was possible to vecto- 
rially add the measured winds to the Earth-relative velocity of the KAPPA program as 
they were in the CSA-RADFAST program, convert the drag acceleration to drag force, 
and attain another drag-coefficient variation. This variation is shown in figure 10 with 
CSA-RADFAST R4 and reference 1 drag curves. The basic shape of the KAPPA curve 
is similar to the others, but its drag rise occurs at a higher Mach number. 
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Figure 8.- Drag results obtained by the application of the 
CSA-RADFAST method to models of the test series 
compared with the results of reference 1. 
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Figure 10.- R4 drag results obtained by two different radar -data -alone 
methods compared with the results of reference 1. 
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ANALYSIS OF RESULTS 


In attempting to give a physical explanation for what caused the errors in CSA- 
RADFAST results, three possible error sources may be immediately identified: (1) the 
data reduction equations, (2) the smoothing techniques, and (3) the radar measurements. 

A six-degree-of-freedom trajectory program (ref, 5) with a radar subroutine was used 
to generate radar data. In addition to the four quantities (p,^ r ,0 r ,p D ) usually obtained 
from actual radar tracks, the radar subroutine also yields >p r and The main pro- 

gram simultaneously computes the total aerodynamic force coefficient Crp the compo- 
nents of which it is actually using internally in computing the trajectory. It should be 
remembered that for small angles of attack, C T « C D < The output of this generated 
case was designated GEN 1. In order to verify the data- reduction equations, the first, 
derivatives (p, ij/ r , and 0 r ) of GEN 1 were numerically differentiated to yield p, ^ r , 
and 0 r . These nine quantities (the three just obtained and the six from GEN 1) were then 
used as input to a RADFAST program modified to accept them directly. This checkout 
case was designated GEN 2 and its output, C T versus M, is shown compared with the 
C T versus M of GEN 1 in figure 11(a). It can be seen that there is extremely good 
agreement. 

The next step taken to help gain an understanding of the problem was to use just the 
four quantities (p, i// r , 0 Y , and p D ) from GEN 1 and to treat them exactly in the 

same manner that actual radar data (which, in fact, is comprised of only these four 
quantities) is treated. In accordance with the procedures described in the section enti- 
tled '’Analytical Methods” these generated radar data were smoothed using the CSA/CRT 
technique. Its output, the coefficients to the cubic equations, were used as input to 
RADFAST where they were evaluated and used to compute a C^, versus M curve. 

This output is designated as GEN 3 and is shown compared with Cy versus M of 
GEN 1 in figure 11(b). It can now be seen that while some degradation has occurred the 
agreement between the curves is still good. 

Actual radar flight data has several things which could affect the accuracy of the 
reduced Cjy the basic calibration of each measurement (bias, linearity, etc.), the 
sampling rate (Is it high enough to define all frequencies of interest?), and noise (white, 
black, signal-to-noise ratio, etc.). For the purpose of this study, only the effect of a 
fixed sampling rate and white noise with a standard deviation comparable to that found on 
the Wallops FPQ-6 radar is examined. First, the effect of the fixed sampling rate alone 
is studied and then the effect of a combination of sampling rate and noise is investigated. 

The six-degree-of-freedom trajectory program has a computing interval which is 
controlled internally by an error-estimation criteria. In effect, the generated radar data 
turns out to have a variable sampling rate. However, since the variable computing 
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(a) Checkout case of the data reduction equations. 



(b) Checkout case of the CSA/CRT technique. 

Figure 11.- Results of a study to identify the factor most responsible for 
' the flight results. 
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(sampling) interval is generally smaller than the fixed sampled interval of the radar, it 
is a simple matter by interpolation to produce a new set of data with equal intervals. 
Again, using the procedures of the section "Analytical Methods" another Crp versus M 
curve, designated GEN 4, is achieved and is compared with that of GEN 1 in figure 11(c). 
Since the agreement is still good, it appears that a fixed data-sampling rate has a small 
effect on the data reduction. 

White noise is random in nature, has a Gaussian (normal) distribution with a mean 
of zero. Quasiwhite noise was generated by a FORTRAN IV subroutine (GETRAN) with a 
standard deviation comparable to that found on actual radar data and algebraically added 
to the generated radar data of GEN 1. The following table presents the estimated stand- 
ard deviation found on the Wallops FPQ-6 radar for these flights: 


p, m . . . 

deg . 
Q v , deg . 
P D , m/sec 


. . 0.9144 

1.123 x 10" 3 
1.123 x 10" 3 
. . 0.2253 


By interpolation, a new set of noisy data with equal intervals was generated. And finally, 
once again using the data-reduction procedures, another C T versus M curve desig- 
nated GEN 5, is computed and compared with the C-p versus M curve of GEN 1 in 
figure 11(d). Now it is possible to see a degradation in the results comparable to that 
attained when using actual radar data. 

Thus, in spite of employing advanced meteorological measurement techniques, 
optimizing the trajectory, and updating both the data reduction equations and the filtering 
procedures, the precise definition of the drag-rise Mach number from radar data alone 
is not possible because of the inadequacy of the cubic spline approximation technique to 
handle the noise inherently found on the Wallops tracking radars. 


CONCLUDING REMARKS 


An investigation was conducted of new techniques used to determine the complete 
transonic drag characteristics of a series of free-flight drop-test models using, for the 
most part, radar -tracking data only. This work was done to improve the results of that 
which is typically a back-up mode to the accelerometer technique used in free-flight drag 
studies. In this effort, exact kinematic equations tailored to enhance smoothing proce- 
dures were formulated. Preflight trajectory design optimized the contribution of the data 
known to be most accurate. The latest meteorological measurement systems were 
employed. Analytical filtering techniques supplemented by visual displays which permit 
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on-the-spot engineering judgment to be exercised in the postflight analysis of the data 
were used. The results of this exploration were compared with those obtained by a more 
direct, accelerometer approach using onboard, instrument data. The comparison indi- 
cated that the accelerometer drag curve was approximated by the radar-data-alone drag 
curve. A small amplitude perturbation on the latter curve, however, precluded a pre- 
cise determination of its drag rise. On the other hand, the radar-data-alone approach 
has merit as a backup method, for the drag rise could be estimated. In an attempt to 
define physically which factor caused the above results, an analytical study using gener- 
ated radar data was made. In this study it was found that the filtering techniques 
employed were not capable of accounting for radar noise. When the generated data were 
degraded by white noise comparable to that inherently found on the Wallops FPQ-6 radar, 
the results were essentially the same as when using real radar data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 11, 1974. 
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APPENDIX A 


HAND -CALCULATION DAT A -REDUCTION METHOD 


The most significant equations of the method vised in the past for determining the 
variation of drag coefficient with Mach number from principally radar data took the form: 

Drag force along the flight path: 



W / dV E 
6 \ dt 


+ 


g sin 


(Ai) 


Earth-relative velocity: 



Free-stream velocity: 

V » * V E + V w cos ( y E,y - * w) cos *E,p (A3 > 

Drag coefficient: 


C - D 

C D- cQ? 


(A4) 


These equations were formulated based on the following assumptions: 

(1) Ballistic flight vehicles would fly at small angles of attack and consequently all 
aerodynamic force would be drag force. 

(2) Other than the drag force, the only force acting along the flight path would be a 
component of the model weight. The acceleration proportional to the algebraic sum of 
the above forces is the scalar magnitude of the velocity-vector rate of change. 

(3) The Earth is considered flat and nonrotating due to the close proximity of the 
flight vehicle to the radar and the shortness of flight duration. 

The first assumption is of utmost importance and is examined in some detail in 

/ 

appendix C. The second assumption forms the basis for equation (Al) which sums the 
forces along the radar flight path (the Earth-relative velocity vector). Drag acts along 
the free-stream velocity vector; however, at the test speeds of interest and the low actual 
wind velocities, these two velocity vectors are almost coincident. The first part of the 
third assumption is quite appropriate. The Earth could be considered flat for flights of 
interest here. The second part of the third assumption, however, is erroneous. While 
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the duration of the flights was short (on the order of 10 sec for the helium-gun flight tests), 
the Coriolis forces which arise from tracking a free-flight body using a fixed radar on a 
rotating Earth act on the flight model the entire time, no matter how short, and conse- 
quently influence the force coefficients which are being sought. For low drag models the 
Coriolis forces can be significant. In the original helium-gun test setup (see fig. Al) the 
models were launched directly away from the Doppler radar set and usually into or with 
the surface winds. The models were launched with low elevation angles and quite often 
reached an altitude no greater than 610 m. In many cases the models were flying more 
or less horizontally in the region of interest. Thus, the free-stream velocity could, for 
the most part, be obtained merely by adding or subtracting, depending on its direction, 
the actual wind speed. In equation form, 

v = V™ ±v (A5) 

00 E w 

The sign would be positive for . a headwind and negative for a tailwind. Since the line of 
sight from the Doppler radar set was, for all practical purposes, tangent to the model 
flight path throughout the test region, it was also possible to consider the Earth -relative 
velocity to be composed entirely of the velocity obtained from the Doppler set (V E = p D ). 
Furthermore, the rate of change of this Doppler velocity was the dV E/ /dt used in 
equation (Al). In some of the helium-gun flights and many of the rocket-propelled 
flights, these special flight conditions were not possible and more complete expressions 
had to be employed. The velocity equation was modified to include the transverse compo- 
nents of motion as well as the radial one (see eq. (A2)), Correction factors were applied 
to equation (A5) to account for wind which were other than headwinds or tailwinds and for 
the cases where models were not flying horizontally (see eq. (A3)). The second assump- 
tion, however, where only quantities acting along the flight path are considered, was 
maintained so that equation (A5) was modified only in such a manner as to consider the 
component of the actual wind parallel to the flight path. An even more complete expres- 
sion would have been to vectorially add the winds to the Earth -relative velocity. 

The free use of equation (A5) has been responsible for the erroneous practice by 
some analysts of using the time-differentiated free-stream velocity dV^/dt in 
equation (Al). This error is not immediately detectable because with a constant wind 
velocity and a near constant flight altitude, dV^/dt is equal to dV E /dt. However, while 
the wind may be constant at a given altitude, it does usually change with altitude. If the 
flight vehicle is changing altitude and if the time-differentiated free-stream velocity is 
used, an acceleration is introduced which is equal to the product of dV^jdh and dh/dt, 
where the former is the change in wind velocity with altitude and the latter is the change 
in altitude of the flight vehicle with time. This can be shown in equation form by differ- 
entiating equation (A5) as 
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(a) View of vertical line-of-sight plane of Doppler radar. 


Remote positional radar 

\ 





D °PP ,er radar . Surface wind 


(b) Overhead view. 

Figure Al.- Schematic of the launch model typically used in the past. 
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dV^ dV P dV 
dt dt dt 

where it can be seen that 

dV dV 

v w _ w dh 

dt dh dt 

Obviously, in the flight-test cases of constant altitude the dV w /dt term vanishes. For 
the other cases dV w /dt does exist but is an acceleration that has nothing to do with the 
external forces acting on the flight vehicle. 


Management and Reduction of Data to Coefficient Form 

The hand-calculation method was aimed directly toward determining the variation 
of drag with Mach number. For instance, it computes only enough trajectory parameters 
to correlate the meteorological quantities with the aerodynamic forces. The method is 
presented in a step-by-step manner as follows: 

(1) Plot the raw radar data (^p,^ r ,6 r ,Pj^ versus time (t), hand fair the plots, and 

read smooth values of same 

(2) Compute altitude (h), horizontal range (R), north-south component of horizontal 

range (x), and the east -west component of horizontal range (y) using smooth 
values of p, ip v , and 0 T 


h = p sin 0 r (A6) 

R - p cos <? r (A7) 

x = R cos if/ r (A8) 

y = R sin \p r (A9) 


(3) Plot h versus R, hand fair the plot, and read slopes from the faired plot to 

obtain the flight-path angle in the vertical plane 

= tan - ! Ah_ (A10) 

%P AR 

(4) Plot x versus y, hand fair the plot, and read slopes from the faired plot to 

obtain the flight-path angle in the horizontal plane ^y^ \ 
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y E v = tan_1 7^ 
">y Ax 


(All) 


(5) Plot smooth values of ip and 9 versus t, read slopes from the respective 

♦ 1 l 


plots to obtain <//. and 9 r 




6 r = 


At 

A 8 
r 

At 


(A 12) 


(6) Plot \p and & r versus t, hand fair the plots, read smooth values of ip r 

and $ , and convert the units of these quantities from deg/sec to rad/sec 

(7) Compute the Earth-relative velocity with respect to the radar (Vgj using 

equation (A2) and the above converted ip and 9 r 

(8) Plot Vg versus t, hand fair the plot, and read slopes off the plot to obtain 

the scalar magnitude of the velocity -vector rate of change ^dVg^dtj 

(9) Plot dVg^dt versus t, hand fair the plot, read smooth values of dVg/dt, 

and compute the drag (D) using equation (Al) 

(10) Using wind speed and wind direction (p w J from the radiosonde data, 

compute free-stream velocity using equation (A3) 

(11) Using density and speed of sound (V c ) from the radiosonde data, compute 

dynamic pressure ^q^ and Mach number (M) 


q*> = f p^l 


(A13) 


v c 


(A14) 


(12) Finally, compute drag coefficient using equation (A4). 
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EQUATIONS OF THE NEW DATA-REDUCTION METHOD (RADFAST) 


This appendix presents the equations used herein to determine the kinematic motions 
of a vehicle in free flight over a rotating oblate spheroid and the variation of drag coeffi- 
cient with Mach number (RADFAST). These equations were formulated in such a manner 
that each radar measurement and its applicable derivatives could be used directly in the 
determination of the flight trajectory and the drag-coefficient variation with Mach number. 
The purpose of this approach was to enable any one set of data to be smoothed independ- 
ently of any other set. Thus, after the contribution of a particularly accurate set of data 
as a component of a total desired quantity had been maximized by preflight trajectory 
design, it would not then be degraded by being mixed with much less accurate data and 
then smoothed. This approach also facilitates the introduction of the Doppler velocity 
Pp in the data reduction procedure. 

In addition to updating the techniques for obtaining drag coefficients from radar data 
alone, this program also updates the trajectory program and is compatible with the six- 
degree-of -freedom trajectory program of reference 2 which is widely used by NASA 
Langley Research Center, NASA Wallops Flight Center, and others. 

The equations are written in a straightforward manner and are presented in an out- 
line format as follows: 

I. Determination of the constants associated with the tracking radar station (fig. Bl(a)). 

A. Geocentric latitude of the sea-level subradar point ^^sl) * s determined by 


^c ( or sl) 


tan 



tan 0 G ( O r sl ) 


(Bl) 


where the geodetic latitude of the radar, $ G (ORg^J = <£q(OR), and the polar and equa- 
torial radii of the earth b ahd a, respectively, are known. 

B. The geocentric radius vector of the sea-level subradar point R^ORg^) is 
found by 



ab 

\J^ 2 cob^ c (OR si ) + a2 S in2^ c (OR SL ) 


(B2) 


C. The geodetic radius vector of the sea-level subradar point R^ ^O r S jJ *s given 

by 
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Z I’ Z GN’ Z RE OR 



(a) Meridian plane view showing detailed geocentric and geodetic latitude 
relationships of the radar station OR. 

Figure Bl.- Radar -station relationships with respect to Earth. 
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R G 




sin 0 C (OR gL ) 
sin 4 > g ( 0R gL ) 


D. The geodetic x'adius vector of the radar itself Rq(OR) 
the height of the radar above sea level h R as follows: 


(B3) 

is found merely by adding 


R g (OR) - ( OR sl) + h R 


(B4) 


E. The geocentric radius vector of the radar itself R^.(OR) is determined by 

R c (OR) = Jr c ( OR Sl ) 2 + h l + 2 [ R C ( OR Sl)JM cos * (I 

where 

0 = 0Q ( OR Sl) “ ( OR Sl) 

F. Finally, the geocentric latitude of the radar itself <£ G (OR) is calculated by 


<£ g (QR) = sin 


Rp, (OR) 

Oort sln *g< or > 


(B6) 


II. Determination of the positional displacement of the flight vehicle with respect to 
various axis systems (figs. Bl(b) and Bl(c)). 

A. The position of the flight -vehicle origin relative to the geodetic axis system of 
the radar X 0R (0) G is determined by 


x or(°)g ~ 


p cos cos 4/ v 
p cos 0 r sin i^ r 
-p sin 0 r 


(B7) 


The X OR -axis is directed toward the North Pole and the Z GR -axis is directed 
downward. The symbols Xq R , yQ R) and Zq R are components of x qR^^G‘ 

B. The position of the flight-vehicle origin relative to a rotating, Earth -centered, 
radar axis system X RR (0)£ is calculated by 


X RE^E = Pl] X OR^G + 


R c (OR)cos ^ c (OR) 
0 

R c (OR)sin <t> c (OR) 


(B8) 
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-*r 
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X 


OR 


2 


OR 


(b) Positional displacement of the flight vehicle O with respect to the 
geodetic axis system of the radar station OR. 



(c) Overall view showing the relationship of the radar axis system to the 

Earth-centered axis systems. 


Figure Bl.- Concluded. 
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where |t{J is a simple, time -invariant transformation matrix defined in equation (B48). 
The X„ tr -axis lies in the equatorial plane and passes through the meridian of the radar. 
The Z RE -axis passes through the North Pole. The symbols x RE , y RE , and z RE are 
components of X RE^E* 

C. The position of the flight -vehicle origin relative to a rotating, Earth-centered, 
Greenwich axis system X G ^(0) E is found by 

X GN^E “ P'2] X RE^E 


where JT 2 ] is another simple, time -invariant transformation matrix in equation (B49). 
The Xg N -axis lies in the equatorial plane and passes through the meridian of Greenwich, 
England. The Z GN -axis passes through the North Pole. The symbols x GN , y GN , and 
Zqn are components of X GN (0) E . 

D. Finally, the position of the flight -vehicle origin relative to a nonrotating, Earth- 
centered, inertial axis system Xj(0)j is determined by 

Xi(°)i = [T 3 (tj]x GN (0) E (BIO) 

where £rg(t)J is a time -dependent transformation matrix defined in equation (B50). The 
Zj-axis also passes through the North Pole. The Xj-axis lies in the equatorial plane and 
is related to X R by a phase angle A R . The symbols x R y R and Zj are components 
of Xj(0) r 

Finally, A R is determined by 

A E = A Ej O + W E (* " *0) (B11 ^ 


where o> E is the Earth's angular velocity and t is time. 

III. Determination of the velocity of the flight vehicle with respect to various axis 
systems. 


A. Velocity of the flight -vehicle origin relative to the X GR (0) G 
determined by 


-axis system can be 


V 0R (0) G - x 0R (0) G = [M(D] 



(B12) 


M 

where j_M(t)j is a time -dependent matrix defined in equation (B51). The Doppler velocity- 
P E is introduced at this point. 
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B. Velocity of the flight -vehicle origin relative to the X RE (0) E -axis system is 

V RE^E ~ ^RE^E " [jlJ^OR^G ( B13 ) 

C. Velocity of the flight -vehicle origin relative to the X ( -.j^(0) E -axis system is 

V GN^E = ^GN^E = (^J^KE^E ( B14 ) 

D. Velocity of the flight -vehicle origin relative to the Xj(0)j-axis system is 

VjfO)! = Xj(0), = [t 3 (0]x gn (O) e + [TjMjX^O^ (BIS) 

Note that Vj(0)j is generally referred to as the inertial velocity of the body with 
respect to the inertial axis system. 

IV. Determination of the acceleration of the flight vehicle with respect to various axis 
systems. 

A. Acceleration of the flight -vehicle origin relative to the XQ R (0)Q-axis system is 



»-s 

f 


V 

^OR^G ~ ^OR^G " ^OR'°^G “ &*(*)] 

0 r 

+ [M(t)] 

®r 


1 

1 

1 

>_ 


B. Acceleration of the flight -vehicle origin relative to the X RE (0) E -axis system is 

A RE^E = ^RE^E = ^RE^E ~ [ T l]^OR^°^G ( B17 ) 

C. Acceleration of the flight -vehicle origin relative to the X£,^(0) E -axis system is 

Agn'O'e - W°>E =■ X gn (0) e = [t 2 ]x re (0) e (B18) 

D. Acceleration of the flight vehicle origin relative to the X I (0) I -axis system is 

A,(0)! = Vj(o), = Xj(0)j = [t 3 (0]x gn (O) e + 2jT 3 (t)]x GN (0) E + [T 3 (t)]x 0N (O) E (B19) 

V. Determination of the flight -trajectory positions and orientations (fig. B2). ' 
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(a) Meridian plane view showing detailed geocentric and geodetic latitude 
relationships of the flight vehicle O. 


7 7 7 

I' GN* L RE 



(b) Overall view showing the relationship of the flight -vehicle axis system to 

the Earth-centered axis systems. 

Figure B2.- Flight -vehicle relationships with respect to Earth. 
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A. Geocentric latitude of the flight vehicle, is 


= tan 



4> c 


where Xj, yj, and Zj are components of Xj(0)j. 

B. Longitude of the flight vehicle, A^-, = Aq = A, is 

A = tan 1 " A E 

C. Geocentric radius vector of the flight vehicle, R^(O), is 

R c (0) = ^x j + Yj + Zj 

D. Geodetic latitude of the flight vehicle, <p G , is 


(B20) 


(B21) 


(B22) 


<f) r ~ tan 


-1 




sin 4> ( 


3 lbj cos <p c 


(B23) 


where Pg is an auxiliary function computed by starting with /1 q = 1 and iterating until 
fig is determined by the expression 


a 2 - b 2 


n a 


(^n-l) 


R c (0)]Jcos 2 (j) c + (p. n _ 1 ) 2 sin 2 0 c 


(B24) 


E. Geocentric latitude of the sea level, subflight vehicle point, <p G (®sl) ’ * s - 


^c(°Sl) = tan “ X (\ tan *G] 


(B25) 


F. Geocentric radius vector of the sea level, subflight vehicle point, R^, (°sl) ’ is 


R, 


:(°sl) - 


ab 


(B26) 


ft 2 cos 2 0 c (O SL )" + a 2 sin 2 7 c (0 SL ) 

G. Geodetic radius vector of the sea level, subflight vehicle point, r g(Osl)> is 


r g(°sl) _ r c(°sl) 


sin ^ C (0 SL ) 


sin 4> g 

H. Finally, the altitude above sea level h is determined by 


(B27) 
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h= -M°sl) 

when sin <Pq 4 0 and by 
h = R 0 (O) - a 
when sin $ G = 0. 

The geodetic axis system of the flight vehicle X(0) G may now be defined. The 
X G -axis is directed north, the Y G -axis is directed east, and the Z G -axis is directed 
down along the geodetic radius vector of the flight vehicle. 

VI. Determination of flight -trajectory velocities (fig. B3). 

A. Inertial velocity of the flight vehicle with respect to its geodetic axis system is 


(B28) 


(B29) 


Vi<°) G = 


G2 L 


V°)i 


(B30) 


where jT G 2 jj is a time -dependent transformation matrix defined in equation (B52). 

B. Velocity relative to the Earth of the flight vehicle with respect to its geodetic 
axis system is 


v e (°)g - Vj(o) G - 


0, x 

R g (0) cos 4 > ( 

_ o . 

The magnitude of Earth-relative velocity V E is 


(B31) 



V E^G 


(B32) 


The flight -path angles associated with this velocity are 


y E,p - tan 

" V E^kG 

f: e<°>?g 

+ V E (0)? G 

in the vertical plane and 



V E>y - tan - 1 

rv E (°) 1G ' 

L V E (0) iGJ 



(B33) 


(B34) 


in the horizontal plane. 

C. Velocity relative to the free stream of the flight vehicle with respect to its geo- 
detic axis system is 
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v ML 



Figure B3.- Earth -relative and free-stream velocity vector orientation of the flight 
vehicle with respect to its geodetic axis system. 
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v w (o) G = v E (0) G -w E (0) G 


(B35) 


where W E (0) G is the wind velocity relative to Earth with respect to the geodetic axis 
system of the flight vehicle and where W E is the magnitude of W e (0)q. 

The magnitude of the free -stream velocity is 


V«,= \ vJO) G \ 

The flight -path angles associated with this velocity are 

- V »<°>kG 


(B36) 


v = tan" * 

=o,p 


in the vertical plane and 


y = tan - * 

00 ,y 


,V ~ (0 4 + V ” (0) jG 




(B37) 


V oo(0)iG 


(B38) 


in the horizontal plane. 

VII. Gravitational acceleration at the flight-vehicle origin given in various axis systems. 

A. The gravitational acceleration at the flight -vehicle origin in the geocentric axis 
system of the flight vehicle is defined as 


G x (0) c = 


-G n 


R 


(B39) 


where 


and 


g t - A i 


g r = A o 


R c (0) 


2A f 


R r (o) 


'^7 sin2d>^ - 3j\sin 


(B40) 


r* — 

a 

,2 

1 - A 

f a 1 4 r 3 

R c (0) 

1 

R c (0)J L 


sin2d)^ - 1 


‘] 


+ Ar 


Rr(O) 


y ^35 sin - 30 sin^c^ + 3^ 


(B41) 
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and Aq,A 1 } A 2 are oblate Earth gravity constants and 
Aq = 9.798141 m/sec2 
Aj = 0.016050 m/sec 2 
A 2 = 0.0000148 m/sec2 

B. The gravitational acceleration at the flight-vehicle origin in the inertial axis 
system is 

Gj(0)i = JT C 2 ^Gj(0) c (B42) 


where 



is a time -dependent transformation matrix defined in equation (B53) . 


VIII. Determination of drag coefficient versus Mach number M. 

A. The total aerodynamic force F^(0)j is calculated using Newton's law as 
follows: 


( 0 )j = m 


Aj(0)j - Gj(O) 


- 1 

I 


(B43) 


where m is the mass of the flight vehicle, 

B. Introducing measured values of free-stream density and free-stream speed 

of sound V c , the total aerodynamic force coefficient C-j. is 



(B44) 


where S is the reference area and the dynamic pressure is 


q 


OO 



and Mach number M is 


(B45) 



(B46) 


C. Finally, assuming very small angles of attack during the data period (77 ~ 0), 



(B47) 
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IX. Matrices. 



, -sin 4> q (OR) 
0 

cos <Pq{ OR) 


0 -cos 0 q(OR) 

1 0 

0 -sin </>q{ OR) 



cos A r 
sin Ag 


0 


-sin A^ 0 

COS Aj^ 0 

0 1 



cos A E 


sin A, 



-sin Ag 0 

cos Ag 0 

0 1 


" y OR 



x OR 



X OR z OR 
'/ x OR + y OR 

y QR Z OR 

f~2 T~ 
f x OR + y OR 

" ^OR + y OR 


fOR 

P 

y OR 

P 



(B48) 


(B49) 


(B50) 


(B51) 


where x oR ,y OR’ z OR are com P° nents of X OR^G and 


P * 



+ z 


2 

OR - 


T 




+ A) 
+ A) 


-sin^Ag + A^ 

-COS <pQ cos ^Ag + A^ 

cos^Ag + A^ 

-cos (pQ sin ^A.g + A^ 

0 

-sin (p G 


(B52) 
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Similarly, 


-sin 4> q cos(A e 
[ T C 2l] = ' Sin *C Sln ( A E 

COS 

The transpose of either j 
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+ x ) 

-sin^A E + x) 

-cos <p£ cos^Aj, + x) 

+x ) 

cos^Ag + x) 

-cos <p c sin^A E + X) 


0 

-sin <p c 

Tea!] 

or jT c2 J is denoted by a tilde. 


\ 

(B53) 
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STUDIES TO DETERMINE ANGLE -OF -ATTACK EFFECTS 

One assumption used in the past concerning angles of attack is still applicable to the 
updated techniques presented herein. This assumption says that a stable axisymmetric 
flight vehicle will fly at small angles of attack, and consequently, all aerodynamic force 
will be drag force. A preflight study was made to determine if this assumption was 
valid. Three areas where angle -of -attack effects could occur were investigated: (1) in the 
presence of gusts and winds, (2) following tipoff from the launch aircraft, and (3) encoun- 
tering pitch- roll resonance phenomena. 

In order to evaluate the effects of gusts and winds on the angle of attack, reference 7 
was studied to understand the nature of winds. This reference stated that outside of the 
Earth's atmospheric boundary layer (above approximately 1830 m), in the absence of any 
frontal or local storm activity, and within time periods of up to as long as an hour, winds 
are steady and uniform and have no vertical components. However, there are wind 
shears that, when penetrated by aircraft or missiles, appear to an observer in the air- 
craft to be gusts as one would think of them on the surface of the Earth. But, in fact, 
there would be virtually no gusts in the flight regime of interest. In order to study the 
response of the flight vehicle to the penetration of wind shears, a six -degree-of -freedom 
computer simulation was conducted. Worst-probable (3 ct) wind-shear gradients 

^ = 0.05 were numerically superimposed on a statistical mean wind for 

Wallops Flight Center. (See fig. Cl.) Also shown are vehicle altitude and angle-of- 
attack response time histories. These indicate that the vehicle weathercocks into the 
free-stream velocity vector and reaches a total angle of attack no greater than 0,45 . 
Therefore, the effects of "gusts" and wind upon the accuracy of the Mach number and 
drag coefficient are small. Since the six-degree-of -freedom simulation shows that the 
model reaches a small total angle of attack (0.45°) while penetrating 3cr wind shears, the 
assumption that the model will fly a zero-lift trajectory is reinforced. 

Again using a six-degree-of -freedom program, a short study showed that tipoff of 
as much as 12° initial angle of attack would damp to acceptable values by the time the 
flight vehicle had reached the test region. 

Another study using measured values of fin tip deflection and best engineering 
estimates of center-of -gravity offset were used in six-degree-of-freedom simulations to 
check for possible roll-resonance problems. None were found in the simulations and 
none were seriously expected in flight, for each model was mass balanced. 
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